Quasiparticle spectrum of a type-II superconductor in a high magnetic field with 

randomly pinned vortices 
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We show that gapless superconductivity of a strongly type-II superconductor in a high magnetic 
field prevails in the presence of disorder, suggesting a topological nature. We calculate the density 
of states of the Bogoliubov-de Gennes quasi-particles for a two-dimensional inhomogeneous system 
in both cases of weak and strong disorder. In the limit of very weak disorder, the effect is very 
small and the density of states is not appreciably changed. As the disorder increases, the density 
of states at low energies increases and the ratio of the low energy density of states to its maximum 
increases significantly. 



The interplay between superconductivity and a mag- 
netic field has attracted interest for a long time. For suf- 
ficiently strong fields the Meissner phase is destroyed and 
a mixed state appears in the form of a quantized vortex 
latticed. The superconductor order parameter has zeros 
at the vortex locations, through which the external mag- 
netic field penetrates in the sample. Contrarily to pre- 
vious understanding, the increase of the magnetic field 
intensity, and its associated diamagnetic pair breaking, 
is counteracted at high magnetic fields by the Landau 
level structure of the electrons!. This leads to interest- 
ing properties such as enhancement of the superconduct- 
ing transition temperature at very high magnetic fields 
where the electrons are confined to the lowest Landau 
level! Associated with the zeros of the order parameter 
in real space are gapless points in the magnetic Brillouin 
zona3, which lead to qualitatively different behavior at 
low temperatures and high magnetic fieldso. 

It has been argued that the gapless behavior is re- 
stricted to high fields very close to the upper critical 
line where the so-called diagonal approximation (where 
the coupling between Landau levels is neglected) is valid. 
It has been shown, however, that the presence of off- 
diagonal terms does not destroy this behavioiu and that 
a perturbation scheme on the off-diagonal tersps is pos- 
sible, as long as there are no band-crossingsu. It was 
shown analytically to all orders in the perturbation the- 
ory on the off-diagonal terms that there is always a dis- 
crete set of points which are gapless that are associated 
with coherent propagation of the quasiparticles (so-called 
Eilcnbcrger points). These nodes are associated with 
the center of mass coordinates of the Cooper pairs and 
not to some internal structure like in d-wave supercon- 
ductors. Lowering the magnetic field, a quantum level- 
crossing transition has been found that eventually leads 
to a sapped regime and to states localized in the vortex 
coresErQ. 

On the other hand, the effect of disorder on super- 
conductivity has also attracted interest for a long time. 
In the case of non-magnetic impurities and s-wave pair- 
ing Anderson's theorem states that, at least for low con- 
centrations, they have little effect since the impurities 
are not pair-breakingfl. In d-wave superconductors how- 
ever, non-magnetic impurities cause a strong pair break- 



ing effectiJ. In the limit of strong scattering it was found 
that the lowest energy quasiparticles become localized 
below the mobility gap, even in a regime where the single- 
electron wave-functions are still extendedBa. This result 
has been confirmed solving the Bogoliubov-de Gennes 
equations -with a finite concentration of non-magnetic 
impuritiesliil. However, allowing for angular dependent 
impurity scattering potentials it has been found that the 
scattering processes close to the-,gap nodes may give 
rise to extended gapless regionstd. The case of mag- 
netic impurities inJhe s-wave case also leads to gapless 
superconductivity^. 

The question we wish to address in this paper is if the 
presence of disorder affects the gapless behavior found in 
the low-T high magnetic field regime of s-wave supercon- 
ductors discussed above. The case of a dirty b-ut homo- 
geneous superconductor was considered beforaij. It was 
assumed that the order parameter is not significantly af- 
fected by the impurities and retains its periodic structure. 
It was found that when the disorder becomes stronger 
than some critical value, a finite density of states appears 
at the Fermi surface. 

In general, since the interactions between the vortices 
are repulsive, a lattice structure is more favorable ener- 
getically. However, if pinning centers are present in the 
system, the higher energies of different configurations of 
the vortices may be offset by the presence of disorder. 
The question then arises if the gapless behavior prevails 
in this more general case. If the magnetic field is very 
high, such that the system is in the quantum limit where 
the electrons are confined to the lowest Landau level, it 
has been shown that for an arbitraryj-epnfiguration of ze- 
ros there is at least one gapless pointJia. In this paper we 
will consider a more general case in which the off-diagonal 
terms are included. 

In the mean-field approximation the Hamiltonian of 
the superconducting system can be diagonalized and the 
energy eigenvalues are the-solutions of the Bogoliubov-de 
Gennes equations (BdG)E3 
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u(r) + A(r)v(f) = Eu(r) 
v(r) + A*(r)u(r) = Ev(r) (1) 
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where A(r) is the order parameter and E the energy. 
If A(r) = 0, the solutions are the Landau eigenfunctions 
which, for a two-dimensional system perpendicular to the 
magnetic field, read in the Landau gauge 



ql 



(2) 



where n is the Landau index, L x is the length of the 
system in the x direction, I is the magnetic length given 
by I — hc/eH and H n is an Hermitc polynomial. The 
energy eigenvalues are those of an harmonic oscillator 
centered at y = —ql 2 
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(3) 



where uj c is the cyclotron frequency. Taking L y to be the 
dimension along y we obtain that —L y /2 < ql 2 < L y /2. 
In the presence of A (r) ^ we can use the Landau basis 
like 



*(0 = y^,u nq 4> nq (f) 

n q 



(4) 



nq 



and obtain the corresponding eigensystem 

mq 

- [n - n c ] < fc + J2 < q ( A £n) * = e*<k (5) 



in q 



where e M = E^/(hui c ), n c is defined by /i = fno c (n c + |) 
and 



A^ = / drr nk (f)^-^ mq (r). 



(6) 



The excitation spectrum is then obtained solving eqs. (5) 
with an appropriate choice for the order parameter. 

In the lattice case, Ahrikosov's solution can be written 
in the Landau gauge asa 



A(rO = A^e-^ 2 e i2 ?-e-(^0 2 



(7) 



The vortex lattice is characterized by unit vectors a = 
(a, 0) and b = (b x , b y ) where b x = 0, b y = a for a square 

lattice and b x — |a, b y = -^a for a triangular lattice. 
This form for the order parameter is valid sufficiently 
close to the upper critical field, since it is entirely con- 
tained in the lowest Landau level of Cooper charge 2e. 
We will not consider contributions to the order parameter 
from higher Landau levels. In this case the self-consistent 
equation for the order parameter reduces to a single rela- 
tion between A and V, the attractive interaction strength 



between the electrons. In the following we will consider 
the square lattice (L x = L y = L), for simplicity. In 
this case the lattice constant a — ly/ir and the zeros are 
located at the points Xi = (i + ^)ly/ir, yj = (j + ^)lyfir. 

An expression for the order parameter has also been 
found for an arbitrary distribution of the zeros which in 
the symmetric gauge can be written astZI 



A(x,y)=AH 



x-Xi .y-yi 

: h I : 
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(8) 



Here, is the number of vortices in the system (num- 
ber of zeros (xi,yi)). The solution of the BdG equations 
is then simply obtained numerically using eq. (6) and 
performing the gauge transformation of eq. (8) to the 
Landau gauge. 

In the lattice case it is more convenient to use a 
resentation in terms of the magnetic wave-functionst 
instead of eq. (4), to take advantage of the translational 
invariance. The generalization to a random configuration 
is however more conveniently done using a real space rep- 
resentation. We consider a finite system and we have to 
take account of finite size effects. In particular, the use 
of eq. (8) for A(r) is very sensitive. We have there- 
fore compared eqs. (7) and (8) for a finite system. (It 



is convenient to define new variables X 



L ' 



Y 



and to use that L = l^JirN^,). The effect of the finite 
size was eliminated (less than 1% difference with respect 
to Abrikosov's solution) using periodic boundary condi- 
tions and calculating A(x, y) using eq. (8) over a set 
of zeros contained in a circle of unit radius (recall that 
— h < X, Y < 7}) and relating A with A to give the same 
amplitude. Having obtained the excitation spectrum and 
the eigenvectors we calculated the density of states using 
the expression 
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where the sum is restricted to e M > 0. 

We considered three cases: i) square lattice, ii) weakly 
disordered lattice and iii) randomly pinned configuration 
(strong disorder) . In the case of weak disorder we consid- 
ered a distribution of zeros Xi = Xf + 6(— | + r)/y/N$, 
where Xf is the regular lattice location, 8 is an adjustable 
amplitude and r is a random number < r < 1 (and sim- 
ilarly for Yi). In the strong disorder case we allow Xi, Yi 
to take any values on the system. We calculated the aver- 
age over disorder directly on the density of states. That 
is, we select randomly one configuration of the zeros and 
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FIG. 1. Density of states (DOS) (in arbitray units) as a 
function of energy for n c = 10 ±2, for A = 0.2 for the cases of 
S — 0.01, 0.1, 0.5, 1, 2 and for the case of strong disorder with 
adjusted chemical potential (see text). Only the lowest band 
is shown. 

obtain the excitation spectrum and p(to). We repeat this 
process many times and then we take the average over 
the resulting expressions for the density of states. This 
is the final result. 

To check the accuracy of the method we started with 
the lattice case previously solvedoQ. The Cooper pairs 
are formed from electrons in the same point in space and 
in an energy interval of the order of the Debye energy, 
which is taken to be typically of the order of 10 — 20% of 
the Fermi energy (therefore we consider off-diagonal cou- 
plings between Landau levels with n c —no < n < n c +ri]j, 
where njj ~ u>d/u c )- 

The excitation spectrum-depends on A and on the di- 
mensionality of the latticeO. It also depends on n c but 
appropriate rescalings yield a somewhat universal behav- 
ior for not too large values of AqIj. For very small A (in 
units of the Landau spacing) a diagonal approximation 
gives good results and a gapless behavior is found at a 
series of points in the magnetic Brillouin zone associated 
with the zeros in the real lattice (Eilenberger points) and 
in other points which increase in number as n c grows. 
These gapless points contribute to a density of states 
that vanishes linearly as u> — > 0. As A grows, off-diagonal 
terms have to be included. Their effect is two-fold since 
they affect both the normal self-energy and the pairing 
self-energy. This leads to a shift in the chemical potential 
to simulate the 3D case. In this case it can be showncj 
that the Eilenberger points remain gapless to all orders 
in the off-diagonal coupling, as long as A is small enough 
that no band-crossings occur. For n c not too large 



FIG. 2. Density of states (DOS) (in arbitray units) as a 
function of energy for n c — 10 ± 2, and /i = 10 for A = 0.2 
for the cases of 8 = 0, 0.01, 0.1, 0.5, 1, 2. Only the lowest band 
is shown. 

and with the inclusion of off-diagonal coupling, a pseudo- 
gap opens in the spectrum since the weight of the gapless 
points is small (but nonzero!). (This is peculiar to a 2d 
system since in this case the Fermi surface only passes 
through a finite number of points, in general. In 3d there 
is always a value of k z such that it is possible to find a 
gapless pointlj). As n c grows, and for values of A such 
that no band-crossings occur, the density of states is more 
and more similar to the diagonal approximation. 

We want to see if the presence of disorder affects this 
gapless behavior and therefore we consider first the less 
favorable case where a small density of states is already 
present in the lattice case. We therefore consider a typ- 
ical case of rix-,= 10, no — 2 and A = 0.2 which was 
studied beforeO. We maintained the value of the chemi- 
cal potential fixed and independent of disorder. 

In the lattice case the results are very similar using ei- 
ther the magnetic Brillouin formulation or the real space 
formulation. Due to the finitencss of the system there 
is somewhat more structure in p(u>) but the qualitative 
features are the same: a pseudogap of order 0.02, the lo- 
cation of the maximum around to ~ 0.07 and the width 
of order 0.1 for the lowest band above the Fermi level. 

Considering now the case of weak disorder we per- 
formed averages over 100 configurations for S = 
0,0.01,0.1,0.5,1,2. Our results are shown in Fig. 1. The 
cases <5 = and 5 = 0.01 are virtually indistinguishable. 
As S increases the effect of disorder is clear. For 6 = 0.1 
the width of the band is almost unaltered but the struc- 
ture is now broadened. For larger values of 5 this effect 
is more pronounced and the width of the band increases. 
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In particular, the density of states for small to increases 
and the relative weight of the gapless modes increases 
with respect to the maximum p(io), which remains ap- 
proximately at the same energy. The case of 6 = 2 is 
already very similar to the case of strong disorder where 
the randomness is maximized. Even though the disorder 
strongly affects the u(f) and v(r) amplitudes we found no 
evidence for localization. The lowest energy eigenvectors 
extend considerably over the whole system even though 
they are strongly inhomogeneous. 

The pseudogap of Fig. 1 found in the 2d case for 
n c = 10 with the chemical potential adjusted disappears 
over a wide raiuie,of the order parameter if the number of 
levels increasea-l'El. Also, as discussed above, in 3d there 
is always a value of k z such that the pseudogap vanishes. 
Keeping the chemical potential at /i — 10 the pseudogap 
also becomes small since, even though the Eilenberger 
points are not gapless in general, as long as the order pa- 
rameter is small, the gaps are very small throughout the 
Brillouin zone. In Fig. 2 we consider the case of n c = 10 
and ll = 10 as a function of disorder. The case S = 
shows a small finite density of states at the Fermi level 
(due to the finite size considered) and a large density of 
states at low energies. As the disorder increases its effect 
is very pronounced. The DOS broadens as before and 
extends from zero energy with a zero energy value that 
increases as the 8 increses. For strong disorder (8 = 2) 
the DOS is much larger than the lattice result. 

These results show that in the presence of disorder the 
gapless behavior does not disappear and is actually en- 
hanced. Our numerical results indicate that there is a 
finite density of states at zero energy particularly in 3d 
or if the number of Landau levels is not too small. They 
also confirm that the gapless behavior has a topological 
nature and is not specific to the periodic vortex lattice. 

The author acknowledges helpful discussions with 
Zlatko Tesanovic and partial support from PRAXIS 



Project /2/2.1/FIS/302/94. 



1 A. A. Abrikosov, Sov. Phys. JETP 5, 1174 (1957). 

2 M. Rasolt and Z. Tesanovic, Rev. Mod. Phys. 64, 709 
(1992). 

3 S. Dukan, A.V. Andreev and Z. Tesanovic, Physica C 183, 
355 (1991). 

4 S. Dukan and Z. Tesanovic, Phys. Rev. B 49, 13017 (1994). 

5 Z. Tesanovic and P.D. Sacramento, Phys. Rev. Lett. 80, 
1521 (1998). 

6 M.R. Norman, A.H. MacDonald and H. Akera, Phys. Rev. 
B 51, 5927 (1995). 

7 C. Caroli, P.G. de Gennes and J. Matricon, Phys. Lett. 9, 
307 (1964). 

8 P.W. Anderson, Phys. Rev. Lett. 3, 325 (1959). 

9 K. Ueda and T.M. Rice, in "Theory of Heavy Fermions 
and Valence Fluctuations", Ed. T. Kasuya and T. Saso 
(Springer, Berlin, 1985). 

10 P.A. Lee, Phys. Rev. Lett. 71, 1887 (1993); Y. Hatsugai 
and P.A. Lee, Phys. Rev. B 48, 4204 (1993). 

11 M. Franz, C. Kallin and A.J. Berlinsky, Phys. Rev. B 54, 
R6897 (1996). 

12 S. Haas, A.V . Balatsky, M. Sigrist and T.M. Rice, |cond 
mat/9703082] 

1,5 A. A. Abrikosov and L.P. Gorkov, Sov. Phys. JETP 12, 
1243 (1961). 

14 S. Dukan and Z. Tesanovic, Phys. Rev. B 56, 838 (1997). 

15 Z. Gedik and Z. Tesanovic, Phys. Rev. B 52, 527 (1995). 



P.G. de Gennes, Superconductivity of Metals and Alloys 
(Addison- Wesley, Reading, MA, 1989). 

17 V.G. Kogan, J. Low Temp. Phys. 20, 103 (1975). 

18 Y.A. Bychkov and E.I. Rashba, Sov. Phys. JETP 58, 1062 
(1983). 



4 



